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We obtain an explicit solution for the stationary state populations of a dissipative Fano model, 
where a discrete excited state is coupled to a continumm set of states; both excited set of states 
are reachable by photo-excitation from the ground state. The dissipative dynamic is described 
by a Liouville equation in Lindblad form and the field intensity can take arbitrary values within 
the model. We show that the continuum states population as a function of laser frequency can 
always be expressed as a Fano profile plus a Lorentzian function with effective parameters whose 
explicit expressions are given in the case of a closed system coupled to a bath as well as for the 
original Fano scattering framework. Although the solution is intricate, it can be elegantly expressed 
as a linear transformation of the kernel of a 4 x 4 matrix which has the meaning of an effective 
Liouvillian. We unveil key notable processes related to the optical nonlinearity and which had 
not been reported to date: electromagnetic induced-transparency, population inversions, power 
narrowing and broadening, as well as an effective reduction of the Fano asymmetry parameter. 


INTRODUCTION 

The original Fano model was introduced by U. Fano in 
1935 [T] and formalized in 1961 [5] to explain the asymme¬ 
try in the absorption or photo-current profile as a func¬ 
tion of laser frequency used to ionize a gaz of Helium-like 
atoms [3]. Previously, similar diffuse absorption bands 
induced by iodine-chloride pre-dissociation have been ob¬ 
served and theoretically addressed [3] . Friederichs 0 de¬ 
veloped the mathematical formalism of perturbation of 
linear operators to describe the essential feature of the 
Fano model: a discrete state coupled to a continuum set 
of states; both sets of states being reachable by photo¬ 
excitation from the ground state. The resulting photo¬ 
current, which is proportional to the population of the 
continuum set of states, as a function of the excitation 
laser frequency is known as the Beutler-Fano or Fano 
profile: 


where q is the ratio of the transition dipole moment of the 
ground-discrete and ground-continuum transitions, and 
e = (wl — We )/7 where hwe is the energy of the discrete 
state relative to the ground state, Wi is the incident ra¬ 
diation field frequency, and h'j = mrV'^ is the linewidth 
of the excited state, induced by its coupling (per unit 
of energy) nV^ to the continuum set of states, n being 
the density of states. The Fano literature is extensive 
and we cannot do justice in this paper to all the contri¬ 
butions since 1935. The interested reader is pointed to 
Refs. [GHS]. 

Two important extensions of the original model have 
been considered: inclusion of incoherent relaxation and 


dephasing processes and high field intensities. The mo¬ 
tivation to include incoherent processes was first to de¬ 
scribe the pressure broadening [3] due to elastic colli¬ 
sions, laser phase fluctuations [10] and spontaneous emis¬ 
sion m m- Nowadays, Fano profiles in nanoscale 
structures are standard [13 mi, for example in plas- 
monic nanostructures [TSHni, quantum dots, decorated 
nanoparticles m and spin filters m- The coherent cou¬ 
pling with the incident light induces large Rabi frequen¬ 
cies which in turn compete with the relaxation rates in 
order to modify the stationary state. The ability to pre¬ 
dict the lineshape, and in general to investigate non lin¬ 
ear optical phenomena in the presence of a continuum 
set of states are the main motivations to consider large 
incident field intensity. A growing number of experimen¬ 
tal and theoretical papers have been appearing on the 
subject of coherent control and ultrafast pulses on Fano 
models [301135] . 

In spite of the ubiquity of the Fano interferences, to the 
best of our knowledge, explicit solutions for high laser in¬ 
tensities and including general relaxation processes have 
not been obtained. Even in Refs. [33] [5S], dealing with 
quantum dots, despite some approximations no analyti¬ 
cal expressions are derived that afford a simple physical 
interpretation of the results. 

In a recent work |26j , we investigated the signatures of 
the Fano interferences in the emitted spectrum of a sys¬ 
tem with a vibrational manifold. We obtained explicit 
expressions of spectroscopic observables like Rayleigh, 
Raman and fluorescence emission but restricted to the 
low intensity held limit where the lowest order of pertur¬ 
bation was enough to describe the laser-matter interac¬ 
tion. In this letter, we focus on the description of the 
non-linear Fano effect on the total population of the con- 
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tinuum excited state, the observable in the original Fano 
model, but for the general case of large intensities of the 
laser field and including dissipation processes. Here, un¬ 
like Ref. |26j . obtaining explicit expressions requires non 
perturbative calculations. We present a method that al¬ 
lows us to obtain such an explicit formulation for large 
field strengths and Markovian baths in an elegant and in¬ 
tuitive framework where the entire solution is formulated 
in terms of a 4 x 4 matrix corresponding to an effective 
Liouvillian acting in the space of the discrete states only. 
We calculate the usual Fano observable, that is, the to¬ 
tal population of the continuum set of states, which is 
related to the optical absorption or to the photo-current, 
as a function of the incident laser frequency. The striking 
result is that such a function can be written exactly as 
a linear combination of a Fano profile and a Lorentzian 
function like in Ref. [55], but where the Fano q and e co¬ 
efficients become effective parameters that are functions 
of the field intensity and the decay rates. 


MODEL 

Although the objectives and results of the present work 
differ from our previous one |26| . the model is similar. As 
there are some differences and for the sake of introducing 
the notation in a self-contained way, we summarize it 
below. The ingredients of the model are schematically 
presented in Fig.j^ The Hamiltonian H — Hq+Hv + Hp 
is exactly the same as in the original Fano model |5]: 

Ho = Eo\g){g\ + Ee\e){e\ + J dkek\k){k\ 

Hv = j dkV{k)\e){k\+V^[k)\k){e\ 

Hp = E[ge cos{ujLt)\g){e\ + g-l cos{u}Lt)\e){g\] 

+ F J dk[fj.kCOs{ujLt)\g){k\+fj.lcos{ujLt)\k){gW, ( 2 ) 

where Hq is the site Hamiltonian, Hy is the coupling of 
the excited state to the continuum. For simplicity, in the 
following, we will consider that V{k) = {e\Hv\k) is real. 
Hp is the interaction with the incident radiation field, 
allowing transitions from the ground state to the discrete 
excited state 5 O e and to the continuum of states 5 -f-)- fc, 
g-ij = {i\g\j) is the transition dipole moment between 
states i and j and E is the field amplitude. 

The relaxation and dephasing processes are introduced 
in an analogue way as in Ref. [55]. It consists in a 
Liouville equation in Lindblad form to insure complete 
positivity of the density matrix describing the quan¬ 
tum system. The dynamics of the system is given by: 
^ = C{t)p, where C{t) = Cpit) + , with hCp = 

—i{l H{t) — H{t) 0 1) is the Hamiltonian conser¬ 
vative part, with H the complex conjugate of i/, and 



FIG. 1: Energy levels and transitions of a Fano-type model 
with dissipation (see main text for notations) 

= T^op + ^^ure is the generator of dissipative dynam¬ 
ics. 

L^„p = J dkrik){Aikg)(^Aik,g) 

- ^ [^<»A^{k,g)A{k,g)+A\k,g)A{k,g)<Sil] | 

+ Te^A{e,g)®A{e,g) 

- i [1 0 Ai (e, g)A{e, g) -b A^ (e, g)A(e, 3) 0 l] |, 

(3) 

describes the population relaxation from the | k) manifold 
and from the |e) state to the ground state. 

Lpnre = - 7 es[|e)(e| 0 \g){g\ + \g){g\ 0 |e)(e|] 

- J dkjkg[\k){k\ 0 \g){g\ + \g){g\ 0 \k){k\] 

- J dk-fke[\k) {k\ 0 |e)(e| -b |e)(e| 0 \k){k\], (4) 

describes pure dephasing, that is the dynamics of the 
non diagonal matrix elements of p. A(i,j) = |j)(t| are 
the jump operators and T{k) is the population relaxation 
rate from state |fc) to \g) as is Fe for the |e) popula¬ 
tion. jij is the pure dephasing rate for the ij coher¬ 
ence. As in Ref. [55], we have used the correspondance: 
\l){m\ <->■ |Z) 0 \m) = \\lm) [57]. We use the rotating- 
wave approximation (RWA) on L = and 

remove non-resonant terms such that L can be consid¬ 
ered time independent. [Ip is a diagonal matrix whose 
matrix elements are equal to Eojl for excited(ground)- 
ground(excited) coherences, and zero elsewhere. 

There are a number of phenomena that occur at strong 
fields which are not described by the present model (ATI, 
continuum-continuum transitions). This is in line with 
previous literature where these processes are also ne¬ 
glected from the model. Within these assumptions, the 
only restriction of the model is that the Rabi frequency 
should be much smaller than the laser frequency, which 
is around the two-level system (TLS) transition as we 
consider near resonant processes. This is once more an 
approximation done throughout the literature. This is 
not restrictive since all the high-held effects described by 
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the model appear at field intensities that are a few per¬ 
cent of the intensities at which the RWA breaks down for 
a TLS in the visible and a linewidth of around 0.1 eV. 
In this context, it is the comparison between the laser 
intensity or more precisely the Rabi frequency and relax¬ 
ation rates that determines the limit between the linear 
and nonlinear regime. 


FESHBACH PARTITIONING AND EFFECTIVE 
LIOUVILLIAN 

The continuum states population, f dkpkk, where p 
is the full steady state density matrix can be obtained 
by finding the kernel of the time-independent operator 
— L), that is (Ol ~ L)p = 0. To solve this equa¬ 
tion, we split L in two terms, L = Lq + V where Lq is 
diagonal and V is purely non diagonal, and proceed to 
Feshbach partitioning [28]. For that, we introduce the 
projectors P = \g){g\ -I- |e)(e| and Q = J dk\k){k\, with 
P + Q — 1. The corresponding projectors for the discrete 
and continuum parts in Liouville space are given by: 

T’ = P(g)P; Q = P(g)Q-pQ(g)P-tQ(g)(3. (5) 

This allows us to rewrite the kernel equation as (£1^ — 
L){V + Q)p = 0. By projecting on both V and Q and 
after some algebra (see Appendix A), we obtain 

- L,s)Vp = 0, ( 6 ) 

with an effective Liouville operator, 

= VLV + VVQGqQVV, (7) 

and Qq = (Q(£li — L)Q)~^. Qq, which is related to 
the resolvent of QLQ, is not straightforward to calculate 
unless QLQ is diagonal. To achieve the calculation of 
0Q, we proceed to a sub-partition of the Q subspace un¬ 
til the projected Liouvillian be diagonal (see Appendix 
A). P^lP — Leff acts on the P space only, but its ker¬ 
nel is equal to the projection Pp of the exact stationary 
density matrix, on the discrete subspace. Once the den¬ 
sity matrix Pp on the discrete space has been obtained, 
the density matrix Qp in the continuum subspace can be 
computed through the following equation: 

Qp=QGqQVPp. (8) 

Leff can be thought as a 4 x 4 matrix when Pp is consid¬ 
ered as a column vector with 4 elements. To obtain an 
explicit expression for Leg, the usual wide-band approxi¬ 
mation is employed. It is also in this same approximation 
that an explicit expression was obtained in the original 
Fano problem [2]. It consists in assuming that the pa¬ 
rameters of the model do not depend upon k. From now 
on, we consider this approximation and define Fc = r(fc), 
/ic = p-k as constants. 


After tedious but straightforward calculations (see Ap¬ 
pendix A), the effective Liouvillian Leg defined in Eq. ([^ 
can be written in a surprisingly simple form: 

Leg = -i{t ® Heg - HeS O 1) + ifff, (9) 

where Heg is an effective Hamiltonian and L^g is the 
dissipative part of the effective Liouvillian. We show ex¬ 
plicitly in Appendix B that Leg has a Lindblad form. It 
can thus be considered as the generator of complete pos¬ 
itive evolution. This ensure that its kernel Pp is indeed 
a physical state. Here, we prefer a different presentation 
of the operator that makes it more amenable for com¬ 
parison to the effective Hamiltonian that is calculated in 
the scattering problem where dissipation is ignored. The 
effective Hamiltonian can be written as: 


Heg = PHqP + iLfieid, with (10) 

p 

Hge\d = —{pe - imTVpc){\g){e\ + |e)(5|). 

It contains an Hermitian and an anti-Hermitian part and 
corresponds exactly to the effective Hamiltonian of typ¬ 
ical resonance problems in Hilbert space 1301 • The 
non-Hamiltonian part of the effective Liouvillian is: 


Leg — Lg + Ac + Tpure! with: 


hLe = (2n7rH^ -|- HTe) 


A{e, g) {e, g) 


-^pure 

hL^ 


^ g)A{e, 5 ) 0 1 -k I 0 A'l’(e, g)A{e, g)^ 

-(jeg + mTpl)[\e){e\ ® \g){g\ + \g){g\ (g) |e)(e|] 
2mrVpci\\gg){eg\\ + \\gg){ge\\) 

( 11 ) 


L^ is a dissipation superoperator in Lindblad form that 
describes the population decay with rate 2mrV'^ Ih + Fg 
due to the coupling to the continuum and the natural de¬ 
cay rate Fe, and i® ^ pure dephasing superoperator. 
L^ is an additional part of the relaxation superoperator 
which cannot be put into a Lindblad form. We stress 
that this presentation allows a comparison to the Hilbert 
space solution but that the full Leg operator of Eq. © 
can be put in Lindblad form (see Appendix B). 

Finally, solving Eqs. Q amounts to finding the kernel 
of a 4 X 4 matrix, which can be explicitely done with the 
help of a symbolic calculation software (see Appendix 
C). Then, applying Eq. (1^ along with the normalization 
condition pgg + pee + jdkpkk = 1 gives us all of the 
populations and coherences. 

The results will be given in terms of dimensionless 
quantities and H'y = mrV^ is taken as the unit of energy. 
In addition to the original dimensionless Fano parameters 
e = {ujl — 9 = Pe/n-irVPc, we introduce the 

new parameter H = pcF/2V, which corresponds 

to a dimensionless R^i frequency. Also, all relaxation 
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rates will be given in units of 7 , this amounts to per¬ 
form the following replacement: Fc —)■ rc/ 7 , Fg —> Fe/ 7 , 

leg leg/l- 


indicates the relative weight of the Lorentzian term in 
comparison to the Fano profile. The transformation into 
the above form involves a rescaling of the parameters 
achieved in the accompanying software (see Appendix C). 


GENERALIZED FANO PROFILE WITH 
DISSIPATION 


The main result of the paper is that the population 
of the excited state ric = f dkpkk can always be brought 
back to a Fano profile /, plus a Lorentzian term: 


^c(Ceffj ^eff) — C 


f (Ceff5 ^eff) 


D 


'^eff ■ 


( 12 ) 


where the dependence upon wl is solely contained in 
Ceff = ■ Wefi, 7 eff and qes are effective Fano pa¬ 

rameters that depend on all the parameters of the model 
but ujl- C is a measure of the total population and D 


Simple explicit expressions for the effective Fano pa¬ 
rameters can be given when the relaxation and the de¬ 
phasing rates concerning the |e) state can be neglected, 
that is when Fe = 0 and leg = 0 (see Eq. Q). This 
is often a very good approximation in the context of 
semiconductor quantum dots or in hybrids consisting 
of an organic molecule adsorbed on metallic or semi¬ 
conductor nanoparticles (l/^e ~ nanoseconds, 

h/n'KV'^ « 10 femtoseconds). In that case, the only re¬ 
laxation process consists in the continuum states popu¬ 
lation decay to the ground state and the profile is given 
by a pure Fano profile and the Lorentzian term is absent, 
D = 0 and C = have obtained: 


7eff r, [1 + (g2 + l)L! 2 (f^ 2 [( 2 f ]2 + + 2/F, + 2)] 

7 2112 Ye 

^efF _ q 2 ^ \ . Qeff _ Ff, 1 

7 7 V 2Vfi + Tc)' q 2112 +rc 7 eff’ 


(13) 


We now discuss each one of the parameters as a func¬ 
tion of the Rabi frequency and illustrate them in Figurej^ 
We note that these parameters have a non-linear depen¬ 
dence on the Rabbi frequency 11, or in other words, a non¬ 
linear dependence on the strength of the field. The pref¬ 
actor C which is proportional to the intensity of the field 
for weak fields (linear regime) saturates to C = 1 when 
211^ > Fe. In Fig. we show the normalized Fano pro¬ 
files at two intensities of the field (H = 0.001 and H = 0.1 
for 9 = 5). As the field intensity increases, we see changes 
in all of the Fano parameters. The effective width 7 eff 
increases or decreases (power narrowing or power broad¬ 
ening) depending on the value of the relaxation (see Fig. 

and Eq. Q). The effective asymmetry parameter 
Qeff decreases monotonically with q (see Figure and 
Eq. @). As shown in the inset, the population of the 
continuum, even for modest values of H is significant, 
underlying the importance of a theory which can handle 
non-negligible population in the continuum set of states, 
contrary to the approximations in Ref. [25] . The decrease 
of 9 eff can be thought as a consequence of the saturation 
of the discrete excited state population. The energy shift 
h{uieS — uje) of the discrete state |e), that can be seen 
as an AC Stark shift, has an interesting behavior. For 
Fc < 2 , the shift is negative if 0 < H < a /(2 — Fc)/2 and 


it is positive if H > y^(2 — Fc)/2. Therefore, Fc < 2, 
H = \/{2 — Fc)/2 is a null point. On the contrary, for 
large relaxation rates such that Fc > 2 , the shift will be 
positive for all values of the field (see Fig. il). 

Inclusion of population relaxation (Fg ^ 0) or pure de¬ 
phasing processes (leg ^ 0 ) results into heavy expressions 
of the states populations that we provide in the accompa¬ 
nying software (see Appendix C). The main qualitative 
features remain unchanged except for the appearance of 
a Lorentzian term. 


Until now, we have focused on the stationary popu¬ 
lation of the continuum set of states f dkpkk- Another 
quantity that can be measured is the photocurrent, that 
is the total flow of electrons in the continuum. Assuming 
that all the electrons emitted in the continuum \k) states 
are collected by an electrode, in the limit of low bias 
voltage, the current intensity / can be obtained from the 
stationary populations as / = limr^-).o |e| [55] . 

where e is the electron charge. It turns out that / can 
also be brought into the form of a Fano factor and a 
Lorentzian factor as in equation (12|. Explicit expres¬ 
sions for the Fano and Lorentzian parameters can be 
given with both the population relaxation and dephas¬ 
ing included: 
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FIG. 2: Effect of the field on the Fano profiles and param¬ 
eters: (a) Fano profiles for q = 5 for fl = 0.001 and = 0.1 
(b) 7 efr for different values of Fc (c) qes/<l for different values 
of q. Inset show ric as a function of q and 11. Upper end (red) 
corresponds to Tie = 1 and lower end (blue) corresponds to 
Tic = 0 (d) UTeff for different values of Fc and u)e = 10 


7 eg— (l+Fe) ^ [H^Fe (g^ "b Fg -|- 1) + (1 + Fg) 2re "b 1) (Fg -|- 7 eg + 1) + (1 + Tg)^ (Fg -|- 7 eg -|- 1)^] ^ (14) 

^ [^"^rg(g^ -b Fg -b 1) -b n^(l -b Fe)(g^ -b 2Fg -b l)(Fg -b 7eg)] 

Weff/ 

+ U 7 F^(l + re) ^ [rg + Fg-b 7gg(2Fg-b Fg7gg-b 9^ + 2Fg-b 7gg-b 1)] , (15) 

Weffj 


to 


tr 


eff 

7 


7 


qV? 

1 + Fg 


= —• 


C**' = 2112. 


(16) 


These expressions give an exact description of the scat¬ 
tering problem, as formulated by Fano in its original 
work [ 2 ], but extended to large field intensities and dis¬ 
sipation processes. 

An interesting consequence of the nonlinear Fano 
effect, is electromagnetically-induced transparency 
(FIT) [37] . Indeed, in the absence of discrete state 
relaxation and dephasing (Fg = 7 gg = 0), then 7?*“' = 0 
in Eq. (151 and therefore the continuum population 
goes through zero when Ceff = —q^s (see Eq. (HI)- 
The phenomenon of EIT has been characterized before 
[38l|39]. In the standard scheme, EIT is obtained with 
two laser frequencies, where one acts as the control 
radiation that creates the transparency window while 
the second one acts as a probe. In our case, which is 
non-standard, the same frequency acts as a control and 


probe radiation, and the transparency window arises 
from two pathways whose destructive interference point 
is tunable via its intensity. It can be shown that the 
condition Ceff = — ^eff is equivalent to = 1 -b | either 
for the light absorption or the photocurrent intensity. 
This phenomenon is shown in Fig. [^a for the case 
e = 0, q = 15 for different values of Fg, and q = 15, 
and Fig. ib Fg = 0 for different values of e. This 
phenomenon provides an interesting tool for devices as 
well as a means for determining the system parameters. 
For example, irradiating at the discrete level resonance 
(e = 0 ) in weak field and increasing the intensity until 
the induced transparency is found determines the ratio 
of transition dipole moment /Xg to the coupling V such 
that iic/2U = 1/F. In the presence of pure dephasing or 
of relaxation from the discrete state, the zero becomes a 
minimum but its position does not change appreciably 
(see Figure]^ (left)). One should take care not to confuse 

















6 



Q Q 


FIG. 3: Population in the continuum showing electromagnet- 
ically induced transparency (EIT) for a) e = 0, q = 15 with 
different values of Fe and b) for q = 15, Fe = 0 different values 
of e 

the present effect with the trivial zero of the standard 
Fano profile at e = —q (see Eq. 0). In the case 
presented here, the effect is non-linear as the condition 
is Ceff = —qeS and both of these parameters depend 
non-linearly on the field. The condition for the minimum 
is thus obtained by adjusting the intensity of the field. 


CONCLUSION 

We have obtained an explicit formula solving the orig¬ 
inal Fano problem for arbitrary relaxation processes and 
large radiative couplings. The inclusion of more than a 
single excited discrete state is straightforward as long as 
the couplings do not exceed the energy gap between ex¬ 
cited states. Furthermore, our approach serves as a step¬ 
ping stone for descriptions going beyond the wide-band 
approximation, already discussed in the scattering frame¬ 
work. Finally, there is a wide class of systems where the 
present model is directly applicable opening new horizons 
in the analysis of Fano profiles under intense fields, as well 
as in applications and devices that exploit processes such 
as population inversion and electromagnetically induced 
transparency. 
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APPENDIX A: STEADY-STATE DENSITY 
MATRIX 

We outline the procedure used to obtain the steady- 
state density matrix by finding the kernel of — L. This 
process can be separated into two different equations, 
one for the discrete subspace and one for the continuum 
subspace. The density matrix in the discrete subspace 
Vp is the solution to: 


• = f dfc||fce)(/ce|| -I- ||e/c)(efc|| 

• Qs = J dk J dk'\\kk') {kk'\\ 


P 

gg, ge, eg, ee 


Q 

P 2 

gk,kg 



Q 2 

P 3 

ek, ke 



Q 3 

kk' 


FIG. 4: Partition of the Liouvillian superoperator 


[nL-[VLV+ PVQgQQVP)]p = Q (17) 

The density matrix in the continuum subspace Qp can 
then be calculated: 

Qp=QQqQVVp. (18) 

The derivation of these equations is presented be¬ 
low. Solving them requires the calculaton of the 
kernel of a 4 x 4 matrix, which we will do with a 
symbolic calculator and knowledge of the resolvent 
QQqQ = {z + Q{ilL - L)Q)~^\z=o- In what follows 
we will omit z altogether since all of the expressions 
for the density matrix require evaluation of the resol¬ 
vent for z = 0. There are problems arising from this 
inversion because the operator in parenthesis is infinite 
dimensional with a continuous spectrum. However, we 
can partition this equation successively until we arrive 
at a partition where the Liouvillian is diagonal. The 
resolvent of a diagonal matrix is trivial and this will 
be the cornerstone of the rest of the calculation which 
will use exact resummations of perturbative expansions 
given by the Lippman-Schwinger equation. 


The first partition divides the continuum from the dis¬ 
crete states and the remaining partitions divide the con¬ 
tinuum subspace such that P 2 , P 3 and Q 3 are diagonal. 
A Lippman-Schwinger recursion equation applied to the 
last partition Q 3 and P 3 gives Q 2 GQ 2 Q 2 , a Lippman- 
Schwinger equation applied to the second partition gives 
QGqQ and a third and final Lippman-Schwinger equa¬ 
tion yields the resolvent in the entire Liouville space Q. 

Taking the upper partition as an example, we show the 
expressions to calculate the resolvent. We assume that 
QQoQ^{Qin^-L)Q)-^ andPOoV = L)V)-^ 

are known from the previous recursion step. We also 
write the non-diagonal part as PLQ + QLV = V. We 
decompose the Lippman-Schwinger equation G — Go + 
QoVG into its partitions. We insert 1 = 7^ -I- Q on both 
sides of V and project the entire expressions onto V and 
Q. After some rearranging: 

vgv = vg^v -p vgoV{vvQgoQvv)vgv 
Qgv = QgoQvpgp 

pgQ^pgppvQQGoQ 

QgQ = QgoQ + QgoQvvgvvQgoQ 


In the case of a Fano-type model, we use three partition 
spaces which are schematically represented in Figure]^ 

• P = \ \99){g9\ \ + Use) (sell + ||e5)(e5l| + l|ee)(ee|| 

• P 2 = ^ dk\\kg){kg\ \ + \ \gk){gk\\ 


We can explain more explicitly, why equations (19) 
solves the problem: the first equation give VgV by the 
inversion (or an infinite resumation) in P space-we can 
give the explicit equation, and the other equations give 
the others projections in term of VgV. 


















In our aim to obtain the steady-state density matrix, 
we do not need to calculate VQV but only VVQQoQW. 
The steady-state density matrix is obtained by — 
L)p = 0. We can project onto V and Q: 

V{n^-L)Vp + V{n^-L)Qp = Q 
2(0^ - L)Vp + Q(ili -L)Qp = 0 


which can be rearranged to get: 

[P^L - L)V 

- VVQ{Q{nL - L)Q)-^QVV]p = 0 


( 21 ) 


and 


Qp= [(Q(OL-i)Q)-^QViP]p 


( 22 ) 


We recognize that we can group the terms in Equation 

(23) 


(211 into an effective Liouvillian L^g 

L^g = VLV + VVQQqQVV 


where Qq = (2(11^ — T)2)“^, leading to expression (5) 
and (6) of the main text. To illustrate the procedure 
we show the calculation of the resolvent for the subspace 
7^3 + 23 = 22- We start with the first line of Eq. (19) 
which gives the solution for the resolvent in the V 3 V 3 
subspace. We can rewrite it in the form of an infinite 
series VsGVs = Er=o(^3V230o23ViP35o^3)". 

We denote the argument which is exponentiated inside 
the sum as w^. 


W 3 = — J dkV^g{k, k')g{k, e)||A;e)(A:e|| 

— J dkV'^g{k',k)g{e,k)\\ek){ek\\ 

J J dk'V‘^g{k',k)g{e,k)\\k'e){ek\\ 

f dk f dk'V‘^g{k,k')g{k,e)\\ek'){ke\\ 


(24) 


where g{a,b) = [-■j((ili)a& - Ea + Et) + Tat,] \ and 
Tat, is the dissipative term of the \\ah) element of the 
density matrix. We need to take the geometric series 
of term w^. In the wideband approximation, the only 
terms which will contribute to the final result give for 
Vs{Q 2 {nL - L)Q 2 )-^Vz: 


V3{Q2{nL-L)Q2)-^V3 


^ g(fc,e) 

1 -I- mrV^g^k, e) 
9 {e,k) 

1 -I- mrV'^g{e, k) 


||fce)(fce|| 

\\ek){ek\\ 

(25) 


this amounts to a renormalization which introduces an 
effective dissipation term - and thus linewidth - of mrV'^ 


to the ek and ke coherences due to the coupling of the 
discrete excited state with the continuum set of states. 
This term appears in the final effective Liouvillian as 
a Lindblad dissipation from the discrete excited state 
to the ground state with rate 2mrV^/ti. For a more 
detailed description of the integrals and their evaluation 
see Supplemental Information of |2hj . 


APPENDIX B: LINDBLAD FORM OF THE 
EFFECTIVE LIOUVILLIAN 


We now write the effective Liouvillian (Eq. @) in 
Lindblad form. Lindblad |40] and Gorini, Kowassaowki 
and Sudarshan m showed that the most general form of 
a Markovian process was described by a semigroup. This 
ensures complete positivity and trace-preserving proper¬ 
ties of the dynamical map. The Lindblad-GKS form is: 

^ JV "-1 

L = -i[H, P] + 2 E c.. pE*] + [E.P, F*]) (26) 

* j=i 


where Ei’s form a traceless orthonormal set, iL is a self- 
adjoint operator and Cij is a positive definite matrix. We 
choose as the basis the Pauli matrices: 


1 

0 1 

1 

■ 0 i' 

1 

0 

1—1 

II 

1 0 

II 

-i 0 

q 

N 

II 

-1 

1—1 

1 

0 

_ 1 


(27) 


and express the effective Liouvillian in the new basis. 
The result is: 

Leg = Lo+VVQgQiO)QVr 
L,g-Lo = -^[H\]+L^ (28) 

H' — \/2mrV Pc ffy 


where Lq is the original Liouvillian in the V partition, 
E[' is an (Hermitian) addition to the Hamiltonian and 
is a superoperator in Lindblad form (Eq. (26l). The 
strictly dissipative part can be completely determined by 
the coefficients ca: 


(29) 


the positivity of the Cij matrix we diagonalize it. The 
eigenvalues are: 

1 



C3 

iC3 

C2 

Cij — 

-ics 

C3 

-iC2 


C2 

iC2 

Cl 

7 egj C2 

= Q and 

C3 = 




± 


(H -|- 2 re -|- 7 eg -|- 2 ) 

^(H2_2re+7eg-2)2-8H2 


(30) 


Both eigenvalues are positive proving that the matrix is 
positive definite and that the operator can in fact be 
written in Lindblad form. 
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APPENDIX C: ACCOMPANYING SOFTWARE 

To make the results of this paper more accessible, we 
have included all results, equations and plots in an in¬ 
teractive python notebook. The notebook consists of a 
symbolic part where the expressions that include dissipa¬ 
tion from the discrete state and pure decoherence - too 


intricate to write in an article - can be obtained. The 
second part consists of numerical simulations where the 
plots produced in the paper can be redrawn and where 
any simulation of the solution to the Fano model can 
be plotted. The software can be downloaded from fi- 
noqs.wordpress.com. 



